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Abstract 

The possibility is examined that the masses of the hght quarks u 
and d have varied over the course of the universe's evolution. Such 
a variation would have an effect on axion cosmology, and can be the 
basis of a solution of the so-called axion energy problem. If this is the 
case, a scalar force somewhat weaker than gravity could exist that has 
a range of order 1 km. 

There is a long history of speculation about fundamental "constants" of 
nature varying over cosmic time scales. It began with a suggestion of Dirac 
that time-varying constants might explain one of the large number coinci- 
dences in cosmology [1]. A motivation for more recent speculations is that 
in various kinds of theories the gauge couplings, Gn, and other parameters 
can depend on upon the size of extra space dimensions, or upon the value 
of dilaton, moduli, or other fields. Such theories include Kaluza-Klein theo- 
ries [2], Brans-Dicke [3] and extended inflation models [4], "induced gravity" 
models [5], quintessence models [6], superstring theory and brane scenarios. 
From the empirical side, there have been both claims from time to time of 
evidence of changing "constants" and stringent bounds on such changes. (For 
a thorough review see [7].) 

In order to have "constants" change over cosmic time scales, there typi- 
cally needs to be a field with a very flat potential. Specifically, if such a field 
changes by a factor of order unity over the age of the universe, it should have 
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a mass of order the Hubble parameter, which today is about lO^"^^ GeV. In 
this paper, we look at the possibility that the masses of one or both of the 
light quarks u and d is time-varying. (Most of our discussion will deal with 
the case where both rriu and rrid vary with time. However, most of the effects 
we shall discuss would be the same if only one of them did.) We assume that 
they depend on a field that has a very flat potential, and whose couplings 
are proportional to Mp^. The flatness of this potential is limited by the very 
fact that (j) couples to u and d, since QCD and quark-loop effects generate 
a mass for that is at least of order h?QQp)/Mpi ~ 10~^° GeV. This imphes 
that (j) ought to evolve rapidly once the temperature of the universe cools 
to 0{Kqcd)- One might think that this prevents interesting consequences 
for physics today. However, we shall see that such a time variation of m„ 
and rud can help to resolve the well-known "axion energy problem" [8], i.e. 
the problem that in typical invisible axion models too much energy can get 
trapped in oscillations of the axion fleld in the early universe, unless is 
less than about 10^^ GeV. (A variety of other solutions to this problem have 
been proposed in the literature [9].) 

The point is that the potential energy of the axion field, which is generated 
by QCD instanton effects, vanishes if any of the quark masses vanishes. Thus, 
if at the time when the axion potential was "turning on" (when T ~ ^qcd) 
the light quark masses were much smaller than they are today, the axion 
potential energy would have been correspondingly suppressed. It is crucial, 
of course, that the light quark masses already had their present values by 
the time of Big Bang nucleosynthesis, as otherwise the primordial Helium 
abundance would come out wrong. On the other hand, to help resolve the 
axion energy problem, the light quark masses must not already have evolved 
to their present values before instanton effects became important at T ~ 
^QCD- That means that it is necessary that the field started to change 
significantly at about the time when T ~ ^qcdi not much earlier or much 
later. But this implies a mass for (/> of order t^Q^jj^jMpi^ which is just the 
order that is expected to arise from its coupling to the light quarks, as we 
have already noted. 

Let us take the Yukawa couplings of the u and d quarks to be of the form 

C^Yuk = Y^uZurHu{(I)/Mp,) + Ydd^dRHdi(t>/Mpe). (1) 

We assume that (p couples to no other fields of the standard model, except 
through u and d quark loops. The phase of (p is the axion field, so (p has 
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a non-zero Peccei-Quinn charge. Let us write = p(x)e^^^^\ The radial 
field p is assumed to be sitting now at the minimum po of its potential. 
The values of m„ and that are observed today we denote by m„o and 
771(10- For general values of p, one then has, from Eq. (1), 777„(p) = niuop/po 
and md{p) = mdop/po- The Lagrangian density for p and 9 can be written 
approximately as 

jr{p, 0) - Ip'id.er - \mye' + \{d,pf - V{p). (2) 

The axion mass depends on temperature as well as on the value of the radial 
field p, so we will write ml{p,T) — Ca{T)ml{p,0), where Ca{T) = 1 for 
T ^ ^QCD and Ca{T) falls off rapidly for T ^ ^qcd [10]. The axion mass 
depends on p both because of the fact that a{x) = p{x)9{x) and because of 
the p dependence of m„ and 777^. The axion mass can be obtained from the 
term in the chiral Lagrangian that also gives the neutral pion mass, namely 

-CK(p)e'(^+'^°/-^-) + 777d(p)e^(^-'^°/^-)] + H.c. (3) 

This leads to the zero-temperature potential 

V{p) = -2C(777„(p)+777d(p))+C(777„(p)(^ + 7r7A)^+777<i(p)(^-7r7/,)2) + ... 

(4) 

After diagonalizing the mass^ matrix of 6 and 7r°//^, one has m^{p) = 

2C{mu{p) +md{p))/f^, and 777^(^,0) = 8Cp-2777„(p)777d(p)/[777„(p) + 777d(p)], 

giving 

^2 ^ ^ ^Ifn muip)md{p) 

p2 {mu{p) + md{p)) (777„o + mdo ) ' 
which goes as p~^ for small p, and 

v,mp) = -Kfi f "';^"' (6) 

V THuo + ruda J 

which goes as p"*"^ for small p. Not only the axion mass but also the potential 
for p generated by chiral symmetry breaking depends on temperature. So 
we will write Vqcd{p^T) = —Cp{T)Vqcd{p,Q), where Cp(T) equals 1 below 
Aqcd and falls off rapidly above it. There are also contributions to V{p) 
coming from u and d quark loops. These go as 

Vl-loop ^ -:^{\YX\Hu\'+\Yd\W)p'^VMj,, = --L^{\rnuip)\'+\mdip)f)A', 



(7) 
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where A is a cutoff of order tlie wealc scale. (Tlie same loops contribute to the 
(mass)^ of the Higgs doublets, so that whatever solves the gauge hierarchy 
problem should cut off the loops at the weak scale.) The expressions given 
in Eqs. (6) and (7) are roughly comparable for p ~ po ~ and are the 
leading contributions coming from the coupling of p to the u and d quarks. 
They both give a mass scale for p that is roughly of order AQ^jjj/Mpf. As 
noted above, this is the kind of mass that is required for the mechanism we 
are discussing for suppressing the cosmological axion energy density. Thus it 
is important that there not be much larger contributions to coming from 
other sources, such as supersymmctry breaking. 

The contributions to V{p) in Eqs. (6) and (7) arc both unbounded below. 
It is therefore necessary that there be other contributions to V{p). We as- 
sume that there are corrections to V{p) coming from Planck-scale physics 
that cause it to have a minimum at po ~ Mp£. In particular, suppose 
that there are corrections to the terms in Eq. (1) of order \(f)\'^/Mp^ as 
follows: YjIIURH^(l)il - a„|0| VM|,^) + ^^^^^^//^^(l - aM^Mj,^)- In this 
case, if a„^d are positive, |m„(p)| and \md{p)\ have maxima and therefore 
the terms in Eqs. (6) and (7), which go respectively as —\mu,d{p)\^ and 
— \'>TT'u,d{p)\^ , have minima. To simplify the discussion, let us assume that 
du = dd = d and that V{p) ~ \mu4\^ near its minimum. Then wc can write 
mq{p) = |mqo(p/po)(l — |p^/Po)' Q = u,d, as this leads to a maximum for 
\mq\ (and thus a minimum for V{p)) at po, and gives mq{po) = mqo. Then 
we may reexpress the axion mass as 

m^(p,0)^6 ^^-^'/'^°^"^-^>"°"^^; . (8) 
Pop ("^no + mdoY 

For p -C Po, where Vqcd dominates over Vi-ioop, one has 

yiP,^) = -lrnlf:(£j. (9) 

We are now in a position to disuss the suppression of cosmological axion 
energy. First, let us see in rough terms what happens. The QCD in- 

stanton effects begin to be important when T ~ Aqcd- At that point, 
the Hubble parameter ~ Aq^jj/Mp^, and from Eq. (8), m\{p,T) = 

Ca{T)mlip,0) ~ CaiT)^ ~ c,(T)(A^e^/M|.,)(po/p). Thus, the axion os- 
dilations commence when Ca(T) pi/ pq ^ 1, pi being the initial value of p. 
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The energy density in the axion oscillations at that point is pfm^{pi, Ti) ~ 
{Pi I Po)^QCD ~ {Pi / Mpi)h.QCD- other words, it is suppressed by a factor 
{pi/Mp^Y compared to what it would be in a conventional invisible axion 
model having fa ~ Mp^. Another thing that happens when T ~ Aqcd is 
that the radial field p begins to roll due to the potential energy shown in Eq. 
(9). As p changes, the axion mass changes, both because a = p9 and because 
of the p dependence of m„ and m^. The number of axions does not change 
significantly during this process, however, since the axion mass is varying adi- 
abatically. Consequently, the number density of axions falls as R^^, where 
R is the cosmic scale factor. The mass of the axion "initially" (i.e. when 
axion oscillations start) is of order H ~ AQ^jy/Mpf, whereas by the time 
p has settled at its minimum the axion mass has reached its present value 
of (mTrA/po), which is also of order A'^gco/Mpe. Thus, {ma)f/{ma)i ~ 1. 
Therefore, over the whole period when p is evolving, the axion energy den- 
sity changes by about the same factor as the axion number density, namely 
by {Ri/Rf Y'. It follows that the initial suppression of the axion energy den- 
sity by {pi/MpiY (compared to conventional axion models with fa ~ Mp^) 
will also be approximately the suppression of the final axion energy density. 
Since a suppression factor of about 10~^ is required, it is sufficient to have 
Pi ~ 10~^Mpe, which is of order the unification scale. 

Now we will be more detailed. The equation for the rolling of the radial 
field is 

p+3Hp^-V'{p). (10) 

The condition for when p starts to roll is that p ~ Hp and p ~ H'^p. This 
gives 

m'pi ~ c,{T)mlf^/po, (11) 

where the initial value of p is pi. Calling the temperature at that time Tp, 

defining Tqcd = Vm^/^, and using H"^ = ^ [^gT^^ /Mp^, where g is the 
effective number of polarizations of radiation, we have 

c,(r,rHvW-Jj(^)'g). (12) 

The axion field starts oscillating at a temperature given by 

m^(A, To) = c„(r„)m^(p„ 0) = = y i^^gA /M^,, (13) 
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or 



Ca{Ta)-\TjTQCDy 



135 (Mpt" 



Po 



muonido 



27r3gr \ Po J \p,J {muo + mdoy 



(14) 



Comparing Eqs. (12) and (14), we see that the axion oscillations begin at 
about the same time that the radial field starts to roll, a time when the 
temperature is slightly higher than Tqcd- (We will find that po/Pi ~ 10^.) 
Initially, that is at T = T^, the energy in axion oscUations is {paxion)i ~ 
\ca{Ta)'w?^{pi, Q)pj9^, where 9 is measured from the minimum of the potential. 
By Eq. (13) this is 
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{Paxion)i ~ ~^g~^'^a 



pi 
M 



Pi 



The baryon number density at that time was 



{nB)i 



s 



9- 
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Assuming adiabatic evolution since, one has {nB/s)i — {ub/ s)^ 
r]B/7M, and from Eqs. (15) and (16), 



(m„)/T„ 7.047r / p. 



PB 



{ma)i nip r]B \Mpe 



el 



(15) 
(16) 

(nB/n-y)„o«;/7.04 
(17) 



where nip is the proton mass. However, we know from Eq. (5) that (m^) / 
Dividing by Eq. (13): 



p'o (m.uo+mdo)'^ ' 



{ma)i 



45 \^/^ ^ 2 



Tqcd \ ^/nluQnldo ( Mpc 

Ta J 



muo + rrido \ Po 
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Prom Eqs. (17), (18), and (14), one has 
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^ ^Jniuonido Tqcd 7.04 



T^g rriuo + mdo nip tjb _ 
Using Ca{T) ~ (T/Tqcd)'^ [10], one then has 
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We require that at the present time ^axion/^B be less than about 10. This 
can be achieved with pi ^ 10^'^ po to lO^^po, which is what we found from 
our earlier, cruder discussion. 

We have been assuming that both and depend on p. It is also 
possible to imagine that couples only to u or only to d. The results would 
be qualitatively the same. The axion mass vanishes if any one quark mass 
goes to zero. And the QCD chiral symmetry breaking effects would still lead 
to a linear potential for p given by Eq. (6) with either md{p) replaced by 
rudQ or m„(p) replaced by niuo- 

A potential problem for the scenario of a rolling p field is the energy that 
may get trapped in the oscillations of that field. For a realistic cosmology it 
is necessary that the p particles be able to decay to other massless or nearly 
massless particles, whose energy can be red-shifted away. There is nothing 
in principle to prevent this. 

Suppose that the decay rate of p excitations into massless particles is 
given by Fp = amp, where a is the appropriate combination of dimensionless 
couplings and powers of An. This decay must happen before nucleosynthesis; 
so that a must be greater than about 10~^. Then, between the time p oscilla- 
tions begin and the time they are damped out by decays, the universe expands 
by a factor ^ (a)^^/^, assuming radiation domination. The density of energy 
in the p field when it begins to oscillate is of order |y(po)|- The density of 

radiation at that time is p^ad = -^f^^^p ■~^f^l^(p)l/Po- Thus, initially, 
~ {Po/Mpey, and consequently after the p oscillations decay to massless 
modes the fraction of energy in those modes is Jlpdec ~ {c()~^^'^{Po/Mpey. In 
order for this not to be in conflict with bounds from Big Bang nucleosynthe- 
sis on massless degrees of freedom, one must therefore assume that po is less 
than about which is between about 10^^ Mp£ and lO^^Mp^. 

A crucial question is whether the radial field p mediates a long-range force 
that would violate the equivalence principle. The range of such a force would 
be m;' = [V"(p)]-'/' ~ Po/^/n^) ~ m^%Z.a. ~ 1 k"^- The couplings of 
p to the u and d quarks are suppressed by 1/po, but since this is larger than 
1/Mpi there is obviously a danger of violating current limits on long-range 
scalar forces. 

This danger is neatly avoided in an interesting class of cases, namely 
where all the light quarks whose masses depend substantially on p depend 
on p in the same way. This would trivially be the case if (f) only coupled to u 
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or only to d. If coupled to both u and rf, it would mean that mu{p)/muo = 
md{p) /mdo = /(p)//(po)- (In the example studied above, this would imply 
that the constants a„ and defined before Eq. (8) are equal.) Since V{p) 
depends monotonically on the light quark masses (see Eqs. (6) and (7)), and 
these in turn are proportional to f{p), it is evident that extremizing V{p) 
with respect to p also cxtrcmizcs /(p), and thus cxtrcmizcs the light quark 
masses themselves. That means that if we then expand mu{p)uu + md{p)dd 
about the minimum p — pq, calling the excitations of the radial field p, one 
finds that uu and dd couple to but not p. In other words there is no fifth 
force. In other words, there is a dangerous fifth force only if both m„ and 
nid depend significantly on p and depend on it in different ways, or if V{p) 
does not depend monotonically on the quark masses. For, in either of these 
cases, minimizing V{p) does not extremize the light quark masses. If there 
is a residual fifth force, its couplings to nucleons would be through the gluon 
field, and thus of equal strength for protons and neutrons [11]. 

In conclusion, we have explored the idea that the light quarks may couple 
to a field with a very flat potential and may as a consequence have varied 
over the history of the universe. We have seen that this opens up certain 
interesting possibilities. These include a relaxing of the cosmological bound 
on the axion decay constant, and residual fifth forces. The question arises, 
as in many scenarios with parameters varying over cosmic times scales, why 
the potential of the slowly varying field is so fiat. There is also the question 
what sets the initial value of this field. To relax the cosmological bound on 
fa it must start at a value much smaller than its value at the minimum of 
its potential. 

While this paper was being written, we became aware of a recent paper by 
Berkooz, Nir, and Volanski [12], which also considers the idea of time- varying 
Yukawa couplings of quarks. Their idea was that the quark masses could 
have been larger at the time of the electroweak phase transition in order to 
enhance electroweak baryogenesis coming from the CKM phase. The quark 
masses in their scenario are proportional to some singlet field, analogous to 
the field we call (p. This field gets shifted at the time of the electroweak phase 
transition. This is quite similar to the situation we considered, except that 
we have the shift taking place due to the QCD chiral phase transition and 
suppose the quarks to have been lighter at early times rather than heavier. 
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